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^ ; Abstract 

' An electron beam ion trap (EBIT) is used to measure extreme ultraviolet spectra between 10 nm 

r~| . and 25 nm from highly-charged ions of tungsten with an open 3d shell (W XLVIII through W LVI). 

We found that almost all strong lines are due to the forbidden magnetic-dipole (Ml) transitions 
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within 3(i" ground configurations. A total of 37 spectral lines are identified for the first time using 
detailed collisional-radiative (CR) modeling of the EBIT spectra. A new level-merging scheme for 
compactification of rate equations is described. The CR simulations for Maxwellian plasmas show 
that several line ratios involving these Ml lines can be used to reliably diagnose temperature and 
density in hot fusion devices. 
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I. INTRODUCTION 



Spectroscopy of highly- charged ions of high-Z elements is currently the subject of exten- 
sive research. From a theoretical viewpoint, the accurately measured wavelengths, energy 
levels and transition probabilities provide crucial tests for advanced theories of atomic struc- 
ture in a regime where relativistic and quantum-electrodynamic effects become very strong. 
As for applications, since tungsten is currently considered to be a primary candidate for 
the plasma-facing material in the ITER divertor region [l| , the spectra of its ions in a wide 
range of wavelengths are being studied under various conditions. It is not surprising, there- 
fore, that a large number of research papers on the spectra of tungsten ions measured with 



electron beam ion traps (EBIT) 2|-|l3j. tokamaks 14j-|l6|. stellarators [17|] and other high- 



temperature-plasma devices were published over the last decade. A detailed compilation of 



• isl, 



19|. 



the recent results on spectral lines and spectra of W can be found in Refs. 

In the ITER plasma, the tungsten ions will be transported from the relatively cold divertor 
region to the plasma core with temperatures on the order of 20 keV. Although considerable 
efforts are to be spent to minimize radiative power losses due to emission from highly- 
charged ions of W, very useful information for plasma diagnostics can be derived from 
isolated spectral lines. For instance, the electron temperature can be easily found from 
the ratios of strong lines from different ions through the dependence of the ionization balance 
on Te, and the ion temperature can be derived from the line shapes. Determination of the 
electron density rig from spectral lines, however, is not as straightfoward. Most often it 
involves a comparison of allowed and forbidden lines, and thus this technique relies upon 
knowledge of wavelengths and transition probabilities of the involved spectral lines. At high 
densities, when level populations approach the local thermodynamic equilibrium (LTE), or 
Boltzmann, limit, forbidden lines with transition probabilities many orders of magnitude 
smaller than those for allowed electric-dipole (El) transitions are too weak to be observed in 
the spectrum. In low density plasmas, however, the populations of the excited levels which 
decay only via forbidden transitions can be relatively high and therefore result in strong 
intensities. For each forbidden line there exists a transition range of electron densities where 
electron collisions are comparable to the radiative decay rate. It is in this range of densities 
that one may hope to use a particular forbidden line for density diagnostics. 

Currently, more than 80 forbidden lines in tungsten ions, from W^®^ to W^^"*", are known 



from experimental measurements 20| . The high-multipole lines 



'or instance, in x-rays 




rom W ions were observed, 
, vacuum ultraviolet 



24 1 and ultraviolet (UV) [2|, 3[ ranges of spectra in tokamaks and EBITs. The probabilities 



23| . extreme ultraviolet (EUV) 



of forbidden transitions show very strong increase with the ion charge while the collisional 
damping of spectral lines becomes less effective due to a decrease of cross sections. As a 
result, the forbidden lines are more prominent in the spectra of multiply-charged ions. 

Forbidden transitions in highly-charged ions are also a subject of active theoretical re- 
search with emphasis on their use in plasma diagnostics. The visible/UV magnetic-dipole 
(Ml) J = 2 — 3 line in Ti-like ions was analyzed for density diagnostics in hot plasmas since 
the pioneering work of Feldman et al. Continuing this work, Feldman et al. I27I 1 

performed a systematic study of density-sensitive Ml lines in Ti-like ions and in various 



N-shell ions. Recently Jonauskas et al. 



calculated wavelengths and transition probabil- 



ities of Ml lines in Ad"^ configurations of W ions using large-scale configuration-interaction 



methods. Also, Quinet et al. 



29| performed Hartree-Fock calculations of allowed and for- 



bidden transitions in W I-III, addressing in particular the diagnostics of fusion plasmas. An 
extensive calculation of atomic characteristics of eight isoelectronic sequences of tungsten 
ions in a broad range of wavelengths and transitions was recently performed by Safronova 
and Safronova using the relativistic many-body perturbation theory (RMBPT) j3o|. The 
number of theoretical works on spectral lines and transition probabilities in ions of tung- 
sten is too large to cite here, so we refer the reader to the bibliographic databases at the 
National Institute of Standards and Technology (NIST) for an extensive list of publications 



on tungsten 



3l|. 



An example of a density-sensitive line ratio in highly-charged tungsten ions can be pro- 
vided by the ratio of electric-quadrupole (E2) and magnetic-octupole (M3) lines in Ni-hke 



W^^^ 32|]. These two close lines at about 0.793 nm are due to two 3d^^-3d^As parity- 
conserving transitions. They were experimentally resolved only recently jl2l|. although tne 



and KBIT 
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unresolved spectral feature was known for several years from tokamak [21 
measurements. A large difference in transition probabilities for these lines, on the order of 
10^, results in a different response to collisional destruction of level populations. As was 
shown in Ref. |32|, the E2/M3 line ratio in W^^+ can be used for density diagnostics in the 
range of typical values of rig in tokamaks. 

The goal of the present work is to study the magnetic-dipole transitions within the ground 



configurations of tlie 3(i" ions of tungsten using the NIST EBIT. Previously we reported 
several EUV lines in Co-, Ca- and K-like ions |8|, |9| within configurations 3(P and 3d, 
respectively. Here we extend our measurements to include the remaining ions of tungsten 
with open 3d shell. Using detailed collisional-radiative (CR) modeling, we identify the mea- 
sured spectral lines in the EUV range of spectra between 10 nm and 25 nm. In addition, we 
perform CR simulations to explore potential use of the newly identified lines for diagnostics 
of hot fusion plasmas. 

The paper is organized as follows. Section II describes the experiment and measurement 
of the spectra. Details of the CR modeling are presented in Section III. We then discuss the 
identification of the new Ml lines. Section V presents the analysis of the line ratios that can 
be used for density diagnostics in fusion plasmas. Finally, the last Section summarizes our 
conclusions. 



II. EXPERIMENT 



The measurements of spectra from the 3d^ ions of tungsten were performed at the NIST 



EBIT facility 33| using a grazing- incidence EUV spectrometer 3J]. The photons were 
collected by a spherical gold-coated mirror. A 1:1 image of the EBIT plasma column was 
focused onto the spectrometer entrance slit. The mirror center was at an equal distance of 
480 mm from the EBIT axis and the spectrometer entrance slits. The photons were dispersed 
with a reflection flat-field grating with 1200 lines/mm. The grazing incidence angle for both 
the mirror and the grating was 3°. The slit width was kept at 500 resulting a constant 
resolving power of about 350. The EUV spectra were directly recorded by a liquid nitrogen 
cooled back-illuminated charge-coupled device (CCD) that was placed in the focal plane of 
the grating at a distance of 235 mm. The CCD detector has an array of 1340x400 pixels 
(20um X 20/im each). A detailed description of our EUV spectroscopic system can be found 
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34|. 



The spectra were measured in two separate runs, one in 2008 (run A) and another in 
2010 (run B). The nominal electron beam energies in run A were 4500 eV, 4750 eV, 5000 
eV, 5250 eV, 5500 eV, 6000 eV, and 7000 eV, and the observed spectra were in the range 
between 4.5 nm and 19.5 nm. A theoretical analysis of the spectra indicated that some 
additional lines may have longer wavelengths and therefore the second run of measurements 



was initiated. The beam energies for run B were selected to be complementary to those 
for run A, namely, 4665 eV, 4840 eV, 5155 eV, 5355 eV, 5755 eV, and 6500 eV, and the 
observed spectral window was shifted to 8 nm to 26 nm by translating the detector in the 
focal plane. The electron beam current for both runs was 150 mA and the trap depth was 
approximately 220 V. The trap was emptied and reloaded every 11 s. 

The measured spectra of tungsten were calibrated with lines from lighter elements. Ref- 
erence spectra of Ne, Ar, O, and Fe were measured at several energies between 2 keV and 



9 keV for run A. The calibration of spectra for run B was performed with lines 
O, Fe and Kr at beam energies between 1 keV and 16 keV. Both gas injection 



rom N, 



35| and 



metal vapour vacuum arc ion source (MEVVA) [36| systems were utilized in the calibration 
runs. The observed calibration lines were fitted with the statistically-weighted Gaussian 
line profiles. The calibration curve was a fourth-order polynomial fit of the line centers 
(CCD pixel number) to the known wavelengths. The weighting in the fit to the calibration 
curve was based on the quadrature sum of the statistical uncertainty of our observation of 
the calibration line center, the accuracy of the calibration line wavelength, and estimated 
systematic measurement uncertainty. When a wavelength was measured at various beam 
energies, the final wavelength was taken to be the weighted average of the corresponding 
values (with exceptions noted below), while the total error in the final wavelength was taken 
to be the quadrature sum of the total uncertainty from the calibration curve and the reduced 
statistical uncertainty from the average at various energies. The statistical uncertainties in 
the line positions were typically less than 0.001 nm. The final accuracy of the W spectral 
lines was 0.003 nm. 

The measured spectra for tungsten are shown in Figs. [U (beam energies Eb = 4500 eV 
to 5250 eV) and |2] (beam energies Eb = 5355 eV to 7000 eV). The run-A spectra are shown 
in black and the run-B spectra are presented in red (color online only). The spectral region 
in the figures is limited to A= 10 nm to 20 nm since almost all Ml lines from the 3(i" ions 
are within this range. Only one line, A ~ 21.203 nm in the V-like ion, was found above 20 
nm, and therefore we do not show the run B spectra at longer wavelengths. The identified 
transitions in various ions of tungsten are indicated by vertical dashed lines in the plots. The 
measured spectra also contain a few impurity lines from oxygen (e.g., at 15 nm) and xenon, 
which are marked by asterisks. The highest-energy spectrum of Eb = 7000 eV also shows a 
few lines from Ar- and Cl-like ions which have already been identified in our previous work 



III. COLLISIONAL-RADIATIVE MODELING OF EBIT SPECTRA 



Generally, identification of measured spectral lines greatly benefits from applying meth- 
ods that include comparisons of different physical parameters, such as wavelengths and 
intensities. While atomic structure methods for simple ions can calculate wavelengths with 
the accuracy at the level of 0.01% or even better, the simulations for multi-electron ions with 
open shells may not be as precise as needed for unambiguous line identification. Another 
often used technique in EBIT studies is the analysis of the variation of line intensity with 
beam energy. This method, however, would only be of marginal value when neighboring 
ions do not differ much in ionization energy, and therefore it is difficult to uniquely associate 
a line with a specific ionization stage. 

The most reliable analysis of spectral lines can be accomplished with a collisional-radiative 
modeling of EBIT plasmas. For any given set of plasma parameters, such as beam energy 
and density, the CR simulations can produce a detailed synthetic spectrum containing lines 
from a number of ions. A comparison of calculated line positions and line intensities with 
the spectra measured at several energies provides practically unambiguous identification of 
spectral lines. This method was successfully used in our previous publications jlHl in order 
to analyze and identify dozens of spectral lines from highly-charged heavy ions in x-ray and 
EUV regions. 

In this work we implement the non-Maxwellian collisional-radiative code NOMAD js^l 
for the calculation of spectra from tungsten ions in the EBIT. The solution of the steady- 
state rate equation 

A-N = (1) 

provides populations of all relevant atomic states and, consequently, intensities of spectral 
lines. Here is the vector of populations of atomic states included in simulations and A 
is the rate matrix describing physical processes that affect state populations. The detailed 
representation of Eq. ([1]) is: 



j>i j<i 
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FIG. 1. Tungsten spectra between 10 nm and 20 nm for beam energies between 4500 eV and 5250 
eV. The identified transitions are indicated by vertical dashed hnes. The spectra from run A are 
shown in black and the spectra from run B are s^own in red. Asterisks show the strongest impurity 
lines. 
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FIG. 2. Tungsten spectra between 10 nm and 20 nm for beam energies between 5355 eV and 7000 
eV. The identified transitions are indicated by vertical dashed hnes. The spectra from run A are 
shown in black and the spectra from run B are s^own in red. Asterisks show the strongest impurity 
lines. 



-N., ( E (^S + neR%) + E neRT,. + E neRTu + E ^ei?^:;„. + ^..U.rA = (2) 
\j<i j>i km j 

where Nz^i is the population of atomic state j in an ion z, is the radiative transition 
probabihty, Ue is the electron density, -R^fij) -^ffij' R7-i,ki ^^e rate coefficients for 
electron-impact excitation, deexcitation and ionization, respectively, -R^mi is the rate coef- 
ficient for radiative recombination, tiq is the density of the neutrals in the trap, and Rf' 
is the rate coefficient for the charge exchange (CX) between neutrals and W ions. Unlike 
Maxwellian plasmas, dielectronic capture (DC) is normally neglected in EBIT coUisional- 
radiative simulations since this resonant process requires an accurate match of the beam 
energy with the DC energy. 

Charge exchange between the W ions and neutrals in the ion trap can affect the ionization 
distribution. The Kronecker factor Sn in Eq. ([2]) indicates that in our model the CX connects 
only the ground states of adjacent ions. It is worth noting that since the ion charge is so 
high {z ^ 50), the contribution of double CX may be comparable to the single CX. We 
are unaware of any calculations or measurements of single or multiple CX cross sections 
between highly- charged tungsten and neutral atoms or molecules, and therefore the Classical 



Trajectory Monte Carlo cross section scaling 



acx = z- 10-^^ cm^ (3) 

was used in calculations. In fact, the precise value of this parameter is not very important 
as it enters the rate equations as a factor in the product riQVQacx, where Vq is the relative 
velocity between neutrals and tungsten ions. Neither no nor vq are accurately known for our 
experimental conditions, so that the product rioVo was used as the only free parameter in 
CR simulations. 

The NOMAD code solves the rate equations ([1]) using externally calculated basic atomic 
data. For the present work the energy levels, radiative transition probabilities (up to electric 
and magnetic octupoles) and electron-impact collisional cross sections were calculated with 



the relativistic Flexible Atomic Code (FAC), which is described in detail in Ref. 39|]. The 
relativistic atomic structure (including quantum-electrodynamics corrections) and collision 
methods implemented in FAC are well suited for highly- charged ions of heavy elements. 
Our CR model contains 15 ions, from Zn-like W^^"*" to Si-like W^^^ . Since the ions with 
ionization potentials larger than the beam energy Efy have very small populations, only 



6 to 8 ionization stages were kept for each specific simulation. The atomic states for each 
ion contained single and double {An = only) excitations from the ground configuration. 
Single excitations from the 3/ subshells were included up to n = 5 for most of the open-shell 
ions, and up to n = 7 or 8 for closed-shell (Ni-like) ions or ions with one or two electrons 
above closed shells. The double excitations are included only for An = within n=3 shell. 

[7 -11] all singly- and 



In our previous works on high-Z ions with open s and p shells 
doubly-excited states included in the CR modeling were atomic levels, i.e., the fine-structure 
components. Since open 3d"- shells allow many more permitted combinations of angular mo- 
menta, the total number of atomic levels due to single and double excitations from Ss^Sp^Sd"' 
increases drastically and becomes untractable with available computational facilities. While 
a typical number of levels in CR modeling of 4s^4p" and 3s^3p" ions was on the order of 
1000 per ionization stage, the excitations from Ss^Sp^Sd" can generate 10,000 levels or more, 
and thus the total number of levels becomes prohibitively large. 

In order to reduce the size of the rate equations to an acceptable level, the atomic states 
within each 3(i" ion were divided into two groups. The first group was composed of the 
ground configuration levels 3s^3p^3(i" and singly- and doubly-excited levels within the same 
n=3 shell, i.e., 3s^3p^3d"+\ 3s3/3rf"+S 3s^3p^3d''+^ , 3s3p^3d"+'^, and 3/3^"+^. The levels 
in this first group were retained without modification as the fine-structure components. 
The levels in the second group, i.e.. An > 1 excitations 3s'^3p^3d'^~^kl, 3s^3p^3d"'kl and 
3s3p^3d"kl with k > A, were joined into generalized atomic states, which are referred to 
as the "superterms" below. The procedure of level grouping can be exemplified for the 
3d"~^kl configuration. Each of the atomic levels within this configuration can be described 
by the following set of quantum numbers in jj-coupling (FAC level notations are given in this 
coupling scheme): (((3rf" )j_, {3d\)j^)j^, or simply {{j^, jk)j where a+b = n-1, 

j_ and j+ are the momenta of the 3d_ and 3d^ sub-shells, jc is the total angular momentum 
of the core 3d'^~^, = / ± 1/2 is the momentum of the optical electron kl, and J is the total 
angular momentum. Here and below we use notation l± for an / electron with j = / ± 1/2. 
The simplest procedure in level grouping would be to join them according to the atomic 
jj-terms j+)j^, j^). However, the reduction in the total number of states is rather small. 
Even the next level of grouping, based on the core momentum jc, results in several thousands 
of states per ion. Therefore the excited levels in the present work were joined according to 
the j_|_) pairs. For the high excited configurations with a hole in the 3s or 3p subshell, the 

10 



levels were combined by the three momenta {jh,j-,j+) where jh is the hole momentum. We 
found that while such grouping significantly reduces the total number of states per ion, the 
resulting set of superterms provides a sufficiently dense representation of atomic structure 
for each of the 3d^ ions in our CR model. The actual reduction in the number of states in 
an ion can reach an order of magnitude: for instance, for the V-like ion with 3d^ ground 
configuration the originally generated 10801 levels are reduced to 791 fine-structure levels 
and 465 superterms only. This method of including lowest atomic levels and highly-excited 



generalized states is similar to the recently proposed hybrid CR models [40| for high-density 
plasma kinetics, where the highly-excited levels are combined into even more general groups 
of states, namely, configurations or superconfigurations. 

Another feature of our calculations is the additional correction of calculated energies for 
the 3(i" levels. For each of these ions we performed another calculation with the FAC code 
including, in addition to configurations mentioned above, all possible excitations within the 
n=3 complex. (Obviously, the total complex was already included for 3(P and 3d^ ions.) The 
energies of the ground state configurations in the CR model were then replaced by the newly 
calculated energies which, although different by a fraction of a percent only, still improved 
the agreement with the experimental energies and wavelengths. A similar procedure was 



applied in our recent work on EUV spectra from highly- charged ions of Hf, Ta and Au |41 |. 

The transition probabilities and cross sections between the levels and the superterms or 
between the superterms were derived from the FAC results for transitions between atomic 
levels using statistical averaging. The collisional cross sections were then convolved with 
a 45-eV-wide Gaussian electron-energy distribution function of the KBIT beam in order 
to generate the rate coefficients. The final set of rate equations ([T]) was solved in the 
steady-state approximation for a grid of beam energies between 4300 eV and 7000 eV. For 
each energy the charge-exchange parameter UoVq varied between 10^^ and 3 ■ 10^^ cm~^s~^. 
The ionization distributions, level populations, and spectral fine intensities were calculated 
for each combination of Eb and rioVo, and the spectral patterns were compared with the 
measured spectra to find the best agreement. The best value of uqVo was found to be about 
10^'^ cm^^s^^. This agrees with our order-of-magnitude estimates of 10'' cm~'^ for the density 
of neutrals and 10^ cm/s for the relative velocity. 

An example of comparison between the experimental and calculated spectra is presented 
in Fig.[3l The simulated spectrum for the beam energy of 5150 eV with UqVq = 10^^ cm~^s~^ 
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FIG. 3. Comparison of experimental spectrum at the nominal beam energy of 5250 eV (top) and 
calculated spectrum at 5150 eV and nofo=10^^ cm^^s^^. The second order spectrum is shown by 
the shifted line and the strongest impurity lines from Xe (12.4 nm and 13.0 nm) and O (15.0 nm) 
are indicated by asterisks. 



and CX cross section from Eq. ([3]) (bottom) agrees very well with the measured spectral 
pattern at the nominal beam energy of 5250 eV (top); this energy difference is attributed 
to the space charge effects in the trap. The three strong lines at 12.4 nm, 13.0 nm, and 
15.0 nm marked by asterisks are due to xenon and oxygen impurities. Also, Fig. |3] (top) 
shows the second order spectrum shifted along the vertical axis in order to indicate a few 
relatively weak second-order lines. One can see that both line positions and line intensities 
are reproduced in our simulations very accurately so that most of the lines can be identified 
from the visual comparison. Such a good agreement was observed for all cases considered 
in the present work. 
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IV. LINE IDENTIFICATION 



Table H] presents the strongest identified lines between 10 nm and 25 nm in the experimen- 
tal spectra of runs A and B. Almost all lines in this table are the forbidden magnetic-dipole 
transitions within the ground configurations Sc?" of tungsten ions from Co-like W^^"*" to K-like 
W^^"*". The only exception is the 18.468-nm Ml line within the lowest excited configuration 
3p^3(i^ of the K-like ion. Four of the observed lines, namely, 18.567 nm in the Co-like ion, 
17.080 nm and 14.959 nm in the Ca-like ion and 15.962 nm in the K-like ion, are already 
known from our previous measurements |^ . The other lines in Table [T] are reported here for 
the first time. As discussed above, the uncertainty of the measured wavelengths is ±0.003 
nm. 

TABLE I: Identified magnetic dipole lines in the experimental spectra between 
10 nm and 25 nm. The previously known lines are marked by asterisks. The 
FAC level numbers within ions are given in square brackets. Other theoretical 
works: a-Q, b-Q, c-Q. d-Q- 
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[12] ((rf3_)3/2,(4)5/2)4 


3.68(4) 


50 


Cr 


13.886 


13.848 


3d'^ 


[2] 


(4)2 


3d^ 


[15] ((<i3_)3/2,(4)5/2)3 


4.92(5) 


50 


Cr 


14.193 


14.176 


3d<^ 


[2] 


(4)2 


3d^ 


[13] ((<ii)3/2,(4)5/2)2 


3.20(5) 


50 


Cr 


15.363 


15.363 


3d^ 


[1] 


(4)4 


3d^ 


[10] ((rfi)3/2,(4)3/2)3 


1.01(6) 


50 


Cr 


17.133 


17.153 


3d^ 


[1] 


(4)4 


3d^ 


[8] ((rf^)3/2, (4)9/2)5 


6.57(5) 


50 


Cr 


17.826 


17.823 


3d^ 


[3] 


(4)0 


3d^ 


[14] ((rf3_)3/2,(4)5/2)l 


1.32(6) 


50 


Cr 


19.239 


19.317 


3d^ 


[1] 


(4)4 


3d? 


[5] ((rf^)3/2, (4)9/2)4 


3.02(6) 


50 


Cr 


19.684 


19.791 


3d^ 


[1] 


(4)4 


3d^ 


[4] ((rf^)3/2, (4)9/2)3 


2.56(6) 
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TABLE I: Identified magnetic dipole lines in the experimental spectra between 
10 nm and 25 nm. The previously known lines are marked by asterisks. The 
FAC level numbers within ions are given in square brackets. Other theoretical 
works: a-Q, b-Q, c-Q, d-Q. 



Ion 


Sequence X^xp 














A 


charge 




(nm) 


(nm) 


Lower level 


Upper level 


(s-*) 


51 


V 


14.531 


14.511 


3d^ 


[1] 


{d+)5/2 


3d^ 


[9] ({rf'-)3/2, (4)2)7/2 


1.21(5) 


51 


V 


17.215 


17.260 


3d^ 


[1] 


.d+)5/2 


3d5 


[5] ({'^^)3/2, (4)2)3/2 


3.75(6) 


51 


V 


17.660 


17.709 


3^5 


[1] 


d-i- )5/2 


3d5 


[3] ((rf?.)3/2, (4)4)7/2 


1.59(6) 


51 


V 


18.996 


19.098 


3^5 


[4] 


((d3 )3/2,(di)4)ll/2 


3d5 


[13] ((d2 )3,(di)9/2)ll/2 


2.31(6) 


51 


V 


21.203 


21.370 


3d^ 


[1] 


d-)_)5/2 


3d5 


[2] (('il)3/2, (4)4)5/2 


3.40(6) 


52 


Ti 


13.543 


13.521 


3^4 


[5] 


{('i-)3/2,(ti+)5/2)3 


3d* 


[17] ((d2_)o,(d^)4)4 


1.09(6) 


52 


Ti 


16.890 


16.922 


3d-* 


[2] 


(('i?.)3/2>('^+)5/2)l 


3d* 


[7] ((d?. )2,(d2_)4)2 


4.70(6) 


52 


Ti 


17.846 


17.905 


3d-* 


[3] 


{('i?.)3/2i('i+)5/2)4 


3d* 


[10] ((d?.)2,(d2 )4)5 


1.65(6) 


52 


Ti 


19.319 


19.427,19.6'' 


3d* 


[1] 




3d4 


[2] ((dl)3/2,(rf+)5/2)l 


3.31(6) 


52 


Ti 


19.445 


19.568 


3d* 


[3] 


((rfi)3/2,('i+)5/2)4 


3d* 


[8] ((di)2,(d2 )4)4 


3.02(6) 


53 


Sc 


12.312 


12.291 


3d^ 


[1] 


(d''-)-il2 


3d3 


[7] ((d2 )o,(d+)5/2)5/2 


2.75(5) 


53 


Sc 


15.785 


15.812 


3d^ 


[4] 


((d?.)2,(d+)5/2)9/2 


3d3 


[12] ((d_)3/2,(d2 )4)n/2 


1.42(6) 


53 


Sc 


16.027 


16.056 


3d3 


[1] 


m-ii2 


3d3 


[6] ((d2_)2,(d+)5/2)l/2 


1.02(6) 


53 


Sc 


17.216 


17.271 


3d3 


[1] 


i<il)3/2 


3d3 


[3] ((di)2,(d+)5/2)3/2 


2.74(6) 


53 


Sc 


18.867 


18.971 


3d3 


[1] 


(dlh/2 


3d3 


[2] ((d2_)2,(d+)5/2)5/2 


3.41(6) 


54 


Ca 


14.959* 


14.984,15.010'= 


3d2 


[1] 


{dl)2 


3d2 


[4] ((d_)3/2,(d+)5/2)2 


1.81(6), 1.798(6)= 


54 


Ca 


17.080* 


17.147,17.157= 


3d2 


[1] 


(dl)2 


3d2 


[3] (('i-)3/2,(rf+)5/2)3 


3.68(6),3.683(6)= 


54 


Ca 


19.177 


19.281,19.294= 


3d2 


[2] 




3d2 


[6] ((d_)3/2,(d+)5/2)l 


1.72(6), 1.771(6)= 


55 


K 


15.962* 


16.003 


3d [1] (d_)3/2 


3d [2] (d+)5/2 


2. 59(6), 1.48(6)'* 


55 


K 


18.468 


18.536 


3p53d2 [6] ((p^)3/2,(d2 )2)7/2 


3p53d2 [9] (((p3_)3/2,d-)3,d+)9/2 


2.99(6) 



The atomic levels in Table |T] are described in jj-coupling, as calculated by the FAC code. 
The l± groups with total zero angular momentum are not shown. For instance, the excited 
level of the Sc?^ configuration of the Co-like ion has six 3(i+ electrons with momentum 
projections from mj = -5/2 to = -1-5/2 which are omitted in the notation. The numbers 
in square brackets in the level notation columns show the calculated level number within 
the corresponding ion (the ground level is number 1 and so on). 

There are several lines from neighboring ions that have very close wavelengths, e.g., 
18.878 nm in Fe-like and 18.880 nm in Mn-like ion, or 15.368 nm in Mn-like and 15.363 nm 
in Cr-like ion. For such cases the wavelengths were determined from a spectrum where one 
of the lines was strong while the other was weak due to the shifted ionization distribution. 
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The wavelengths for other hnes were obtained by averaging over several measured spectra. 
Table [T] also shows our calculated wavelengths and transition probabilities as well as 



several other theoretical results 



26 



30 



42 



43|. In most cases the present results agree quite 



well with the measured wavelengths although for several lines the difference is rather large, 
as much as 0.5 %. This probably reflects difficulties in atomic structure calculations for such 
complex ions. Most of the calculated transition probabilities are between 10^ s~^ and 5 ■ 10^ 
s~^. The only line with a smaller probability of A = 3.68 ■ 10^ s~^ is the 13.137 nm J=4-. 
transition in Cr-like W. Note also that the recent RMBPT calculations for Ca-like W 
agree with our transition probabilities to within a few percent. 

The energy structure of the 3d"' ions and population flux analysis explain why only the 
forbidden Ml lines are visible between 10 nm {AE ^ 124 eV) and 25 nm {AE ^ 50 eV) under 
EBIT conditions. Normally, the strongest El lines in a coUisionally-dominated spectrum 
are due to transitions between the ground configuration and lowest excited configuration 
of opposite parity. There is, however, a relatively large energy gap between Sp^Sd" and 
3p^3d"'~^^. Our calculations with FAC show that the 3p — 3d excitation energy in tungsten 
ions is about (300-400) eV. Figure H] shows the calculated energy levels below 500 eV for 
Co-like through K-like ions of W. The levels of the ground configuration are represented by 
horizontal lines, and the vertical bars show the spread of the 3p^3d"^^ configuration. For 
W^''^, W'^^"'",W^^"'", and W^^"*", the energy gap between these configurations is larger than 
124 eV, so that the corresponding El lines have wavelengths smaller than 10 nm. For the 
remaining ions, the highest 3d" levels are rather close to the lower edge of 3p^3d"^^ manifold. 
However, only a few possible EUV transitions obey the |AJ| < 1 selection rule for El lines 
and moreover, those transitions are greatly suppressed by small branching ratios due to 
stronger soft x-ray decays into the lowest levels of the ground configuration. 

Some transitions between higher-n states, for instance, n=4-5 transitions in Fe-like and 
Mn-like ions, also fall into the 10-25 nm range. In low-density plasmas, the populations of 
the lowest excited levels of 3d"' are much higher that those of the high-excited states and 
therefore only Ml lines are strong. When the density is high, the populations approach 
the Boltzmann values which are of the same order for all levels in an ion. Since El rates 
between the higher states are much stronger than the Ml probabilities, only the El lines 
will be present at high densities. 

As mentioned above, four of the spectral lines in table [T] have already been observed in 
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Ion charge 

50 51 52 




3d 



3d 



3d 3d 3d 3d 3d^ 

Ground configuration 



3d 



3d 



FIG. 4. Calculated energy levels of the Sp^ScF' configurations in W^''"*" through W^^"*". The vertical 
bars at the top show the spread of the 3p^3(P~^^ configurations below 500 eV. 



our previous experiments. While the new wavelengths for the lines from Ca- and K-like 
ions agree with the known values within experimental uncertainties, the 18.567±0.003 nm 
wavelength for the Ml transition in Co-like ion is shifted with respect to our previous value 
of 18.578±0.002 nm 1^. To address this problem, we reexamined the 4228 eV spectrum of 
Ref. lo| where this line was identified for the first time. As was pointed out in the original 
publication, the line was strongly blended by third-order lines; our current analysis shows 
that its wavelength should therefore have been assigned a larger uncertainty. For the lowest 
beam energies of the present experiment, the Ml line in the Co-like ion is the strongest in 
the spectrum (Fig. [1]) so that its wavelength was determined very reliably. Therefore, the 



presently measured wavelength of 18.567±0.003 nm replaces the previous value of Ref. 



io| 
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Our new wavelength agrees better with the semi-empirical wavelength of 18.541±0.032 nm 
44|. 



V. DIAGNOSTICS WITH THE Ml LINES 

The diagnostic potential of the Ml lines is based on several features. First, the intensity 
ratios for lines from different ions can conveniently serve as a diagnostics of temperature 
and ionization balance. It is also helpful that the spectral window for these lines is rather 
narrow, which reduces dependence of spectrometer efficiency on wavelength. Finally, and 
most importantly, these forbidden lines can be used to diagnose electron density in fusion 
plasmas. 

Various methods have been developed for spectroscopic diagnostics of electron density 



45| . Such techniques make use, for instance, of line intensity ratios or collisional widths of 
isolated lines. The line ratio diagnostics is normally based on comparison of allowed and 
forbidden lines which are populated by similar mechanisms (normally, by excitation from 
the ground state) and whose radiative decay rates differ by orders of magnitude. For low 
densities, when collisional depopulation rates are much smaller than either of the radiative 
rates, the intensities of both allowed and forbidden lines vary linearly with density and 
therefore their ratio is independent of rie. When collisional rates become comparable with 
the probabilities of forbidden transitions, the ratio shows sensitivity to ne, typically over one 
or two orders in ng. A well-known example is the resonance-to-intercombination-line ratio 



46| 



in He-like ions which has been widely used in plasma diagnostics 

While the measured low-density spectra from the 3c/" ions of tungsten contain no strong 
El lines, the decay rates for the observed Ml lines vary by as much as two orders of magni- 
tude, and hence one may expect at least some sensitivity of line ratios to density variations. 
In order to analyze the rze-dependence of the Ml line intensities under typical conditions 
of hot fusion plasmas, we performed another set of calculations with NOMAD using a 
Maxwellian electron energy distribution function and including dielectronic recombination 
within the Burgess-Merts-Cowan-Magee approximation j^. The steady-state solutions of 
the rate equations were determined for electron densities in the range of rif, = 10^*^ — 10^^ 
cm~^ at electron temperature s ~ Iz where is the ionization potential of the ion under 
study. It was recently shown |48| that unlike the low-Z elements, the typical temperatures 
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of the maximal abundance for highly- charged W ions are on the order or even larger than 
the corresponding ionization potentials, and therefore the condition ^ 1^ is well justified 
for hot steady-state plasmas. Also, since the excitation energies of the levels within ground 
configurations are much smaller than T^, the intensity ratios for the EUV lines would only 
be marginally sensitive to electron temperature. Finally, since the intensity ratios involve 
only lines from the same ionization stage, the conclusions will not depend on ionization 
distribution. 

We have already mentioned above that the collisional depopulation of the upper level is 
the primary physical process resulting in density sensitivity of a spectral line. The effect of 
collisions on level population can be approximately parameterized by their fraction in the 



total depopulation rate (see also |46|): 



(^iM = ^ " ^ (4) 



where J2 n^Rlf includes all collisional processes depopulating level such as excitation. 



icol 

deexcitation and ionization. The low-density coronal limit corresponds to — )■ 0, and in 
the high-density Boltzmann (LTE) limit aj — ?• 1. The transitional region for an atomic state 
can be defined by the condition 

0.1 < a.K) < 0.9, (5) 

which covers about two orders of magnitude in Ue, as follows from Eq. (jl]). This region of ne 
determines the range of densities for a spectral line that are most promising for diagnostics 
when compared with other lines that are still in a coronal limit. Figure [5] shows the calculated 
transitional regions of rif. for the upper levels listed in table |T] including, for completeness, the 
levels of Co-like 3d^ and K-like 3d ions. It follows from this plot that, for instance, level 9 in 
the Mn-like ion, from which the 15.368 nm line originates, decays only radiatively {a < 0.1) 
at densities smaller that 2 ■ 10^^ cm~^, while electron collisions are the dominant depletion 
mechanism for this level above 1.5 ■ 10^^ cm~'^. One can see that for almost every ion there 
exist levels with rather different ranges of variations, and therefore one may expect to find 
several pairs of lines whose intensity ratio may be used for density diagnostics. Moreover, 
the transitional regions for a large number of levels overlap with the typical range of in 
core fusion plasmas (marked by the vertical dashed lines). 

The best ne-sensitive line intensity ratios for ions of tungsten from W^^^ to W^^"*" are 
presented in Fig.[6]and are discussed below. The line intensities were defined as J = N-A-AE 
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T 1 1 I I I I I I 1 1 1 I I I I I I 1 1 1 I I I I I I 1 1 1 I I I I I I 1 1 1 I I I I I 




Electron density n (cm ) 



FIG. 5. Transitional (0.1 < < 0.9) regions of densities for the upper levels of Ml lines from 
table [H The numbers are the calculated level numbers in a corresponding ion. The range of typical 
electron densities of core fusion plasmas is shown by vertical dashed lines. 

where is the upper level population (in cm~'^), A is the transition probability (in s~^), 
and AE is the photon energy (in J). 

A. Fe-Iike W'^^+ 

Among the four identified lines in the Fe-like ion, the ratio of lines originating from levels 
4 and 6 offers the best sensitivity to density variations. However, the transitional regions of 
rig for these levels, as follows from Fig. are not too different and therefore the line ratio 
would not change significantly. Indeed, as shown in the top left panel of Fig. El the intensity 
ratio for the spectral lines at 18.988 nm and 15.511 nm varies only within a factor of 2.3 
between 3 ■ 10^^ cm~'^ and 10^^ cm~^. Another complication arises from the overlap between 
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FIG. 6. Density-sensitive line ratios for ions of tungsten from W**^"^ to W' 
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the 18.988 nm lines and the 19.047 nm hne in the next Mn-hke ion W^^^. 



B. Mn-like W^^+ 

Although the seven identified transitions in Mn-like W^^"*" offer several pairs of lines 
of potential usage in density diagnostics (top right panel in Fig. [6]), it would be rather 
challenging for experimentalists to reliably isolate most of these lines from the other lines 
that originate from the neighboring ions. It seems that the 18.670/18.276 ratio that increases 
by a factor of 5 between 10^^ cm~^ and 10^^ cm~^ may actually be the best choice since 
these two lines are well isolated in the measured spectrum. 

C. Cr-like W50+ 

Among all 3(i" ions of tungsten, the Cr-like ion offers the largest number of ne-sensitive 
line pairs. This is due to both the largest number of identified lines and the most separated 
regions of density sensitivity for different lines. The population of level 12, which is responsi- 
ble for the strong and well isolated line at 13.137 nm, starts to deviate from coronal behavior 
at electron densities as low as 9 ■ 10^^ cm~^ (Fig. E]). The strongest radiative transition from 
this level has probability of only 3.68 ■ 10^ s~^ as compared with the typical Ml probabilities 
of 10^ s~^ and larger (Table [T]) and therefore this level becomes collisionally depleted at lower 
electron densities. The line ratios involving the 13.137 nm line, presented in the middle left 
panel of Fig. |6l show the strongest dependence on Ug: for instance, the 19.684/13.137 ratio 
of two strong and isolated lines varies by a factor of 30 between 10^^ cm~^ and 10^^ cm~^. 
The ratios with other lines (middle right panel) vary within smaller limits, from 3 to 6. 

D. V-like W5i+ 

The transitional regions for the levels of the V-like ion are mostly between 2 ■ 10^^ cm~^ 
and 2-10^^ cm~'^. However, electron collisions start to depopulate level 9 at much lower 
densities of about 2 ■ 10^^ cm~'^ and thus the line ratios involving a strong isolated line at 
14.531 nm can be very sensitive to electron density. For instance, the ratio 18.996/14.531 
varies by two orders of magnitude, between 0.1 and 10, for the electron density range of 
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(10^3-10^^) cm-\ 



E. Ti-like W^^+ 

Figure [5] shows that the transitional regions of Ug for the 3(i^ levels in Ti-like ion are very 
close, and therefore the intensity ratios for the observed Ml lines do not exhibit significant 
ng-dependence. We only present a single line ratio 16.890/13.543 (dashed line in the bottom 
left panel of Fig. [H]) which only changes within a factor of 3 within the discussed range of 

F. Sc-like W^^+ 

The two line ratios for the Sc-like ion, 18.867/12.312 and 17.216/12.312, show significant 
variation between 10^^ cm~^ and 10^^ cm~^. The former ratio is not monotonic, as is seen 
from Fig. [6l the low-density limit of about 2 drops to approximately 1.3 before climbing 
to the value of 10 at the high-density limit. This dip results from the presence of other 
metastable levels of 3d^ that interact with the upper levels of those lines. 

VI. CONCLUSIONS 

In this paper we presented measurements and identifications of forbidden magnetic-dipole 
lines within ground configurations of all 3^" ions of tungsten, from Co-like W^^"*" to K-like 
W^^^. Two sets of EUV spectra between 10 nm and 20 nm, independently measured two 
years apart, excellently agreed with each other, thereby confirming a very good reproducibil- 
ity of the results. A total of 37 new spectral lines were identified in the spectra. 

The identification of the observed Ml lines was based on extensive collisional-radiative 
modeling of spectra from non-Maxwellian plasma of EBIT. We introduced a new scheme for 
level grouping in order to reduce the total number of states in our CR model to a tractable 
level. The calculated spectra agree very well with the measured ones, thereby providing 
unambiguous determination of line identifications. 

The importance of the measured magnetic-dipole lines for spectroscopic diagnostics of hot 
plasmas stems from several factors. First, the lines are located within a rather narrow range 
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of wavelengths, which facihtates their measurements and reduces dependence on variation of 
spectrometer efficiency with wavelength. Second, most of the lines are well isolated with only 
a few overlapping to some degree. Third, the intensity ratios of spectral lines from different 
ions can be used to infer electron temperature and ionization balance over a large range of 
plasma density. Finally, as was shown here with a detailed collisional-radiative modeling of 
Maxwellian plasmas, a large number of line intensity ratios are sensitive to electron density 
in the range of magnetic fusion devices. All these features make the Ml lines in Sd" ions of 
tungsten especially useful for plasma diagnostics. 
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